Inelastic electron backscattering in a generic helical edge channel 
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We evaluate the low-temperature conductance of a weakly interacting one-dimensional helical 
liquid without axial spin symmetry. The lack of that symmetry allows for inelastic backscattering 
of a single electron, accompanied by forward-scattering of another. This joint effect of weak interac- 
tions and potential scattering off impurities results in a temperature-dependent deviation from the 
quantized conductance, 8G oc T 4 . In addition, SG is sensitive to the position of the Fermi level. We 
determine numerically the parameters entering our generic model for the Bernevig-Hughes-Zhang 
Hamiltonian of a HgTe/CdTe quantum well in the presence of Rashba spin-orbit coupling. 

PACS numbers: 71.10.Pm, 72.10.Fk 



A key feature of a 2D topological insulator (TI) is the 
presence of gapless edge states at its boundaries with a 
"normal" insulator or the vacuum [TJ [2] . If the system 
is time- reversal (TR) invariant, the counter-propagating 
states of the same energy (carrying momenta k and — k) 
form a Kramers doublet, which makes elastic backscat- 
tering off a potential scatterer impossible [31 0]. Thus, 
potential scatterers on their own cannot prevent electrons 
from ballistic propagation along helical edge states. The 
result is a quantized, temperature-independent universal 
conductance of G = e 2 /h per helical edge. 

On the other hand, inelastic scattering due to a com- 
bination of electron-electron interactions and a potential 
which violates translational invariance along the helical 
edge may affect its conductance; weak interactions lead 
to a temperature-dependent correction reducing the con- 
ductance compared to its universal value. The existing 
theories [3HZ] predict a power-law temperature depen- 
dence (unlike in the quantum Hall effect |5]) and apply 
to interacting helical edges with conserved S z component 
of the electron spin. In the presence of such an auxiliary 
symmetry, the lowest-order processes affecting the con- 
ductance involve backscattering of electron pairs. Such 
two-particle backscattering may result from the presence 
of a lattice potential (Umklapp process) or from an in- 
homogeneity violating the translational invariance of the 
helical edge. In the former case, the temperature de- 
pendence of the leading correction to G is SG oc T 5 , if 
the Fermi level is tuned to the TR invariant point of the 
electron spectrum [3] . In the latter case, the two-particle 
backscattering off an impurity results in SG oc T 6 [6], 
as follows from a straightforward phase space argument. 
If the Fermi level is shifted away from the TR invariant 
point, the Umklapp processes require activation energy, 
and are therefore exponentially suppressed at low tem- 
peratures. For the two-particle backscattering off impu- 
rities, one would expect only a weak sensitivity of SG to 
the position of the Fermi level. 

The prediction of a 2D topological insulator state 
in HgTe/CdTe quantum well heterostructures [TU] 



prompted experiments which indeed found a low- 
temperature conductance close to Go for structures with 
the proper quantum- well thickness |111 112) . The minimal 
model [10] of Bernevig, Hughes, and Zhang (BHZ) is a 
block-diagonal 4x4 matrix Hamiltonian acting in the 
space of four bands originating from two spins and two 
orbital states. This minimal model assumes axial and 
inversion symmetry around the growth axis of the het- 
erostructure (z axis) , which carries over to the edge states 
obtained within the BHZ model. The electron spins in 
these helical edge eigenstates are indeed oriented along 
the ±2 direction at any momentum k. 

The assumed axial symmetry of a HgTe/CdTe het- 
erostructure, even if it exists in the original band- 
structure model, may be lifted by a gate-induced electric 
field in the z direction. This results in a fc-dependent 
Rashba spin-orbit interaction (SOI) [T3j, and the S z com- 
ponent of the electron spin is generally no longer con- 
served [13] . We expect the absence of S z symmetry to 
be a rather generic property of helical edge states, which 
can also be realized in other models. 

In this work, we evaluate the correction SG to the 
universal conductance as a function of temperature and 
Fermi momentum (measured from the TR invariant 
point, k = 0) for a generic helical edge. We will show 
that if the temperature is low and the Fermi momentum 
is away from k = 0, SG is dominated by the combined 
effects of interaction and potential scattering off the dis- 
order potential. In the absence of axial symmetry we 
find SG oc T 4 , which is stronger than the aforementioned 
SG oc T e in the -symmetric case. In addition, SG 
acquires a substantial dependence on the Fermi level, in- 
creasing with its detuning from the TR invariant point, 



see Eq. ( 12 ) below. Moreover, inelastic backscattering of 



a single electron (with energy transfer to another particle- 
hole pair) is possible even without involvement of disor- 
der when one of the participating states is at k = 0. Sim- 
ilar to the two-particle Umklapp process in the axially- 
symmetric model [9], these processes lead to SG oc T 5 if 
the Fermi momentum is tuned to the TR invariant point, 
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and to a thermal activation law for SG(T) if the Fermi 
momentum is tuned away from k = 0, see Eqs. ( 10 ) and 



(111. The crossover between the inelastic backscattering 



involving the k — state and processes utilizing the dis- 
order potential may lead to a non-monotonic dependence 
of SG on the Fermi momentum at fixed temperature. The 
magnitude of the correction to the universal conductance 
and the details of the crossover depend on the specific 
interaction and disorder potentials. However, the very 
existence of the processes we consider rests on the ro- 
tation of the spin quantization axis with k for an ideal 
free-electron helical edge. We determine that rotation 
explicitly by solving numerically the Kane-Mele [9] and 
BHZ 10 j models with added Rashba SOI. The exponents 
of the T-dependence of SG in Eqs. ( 12 ) and ( 10 ) result 



from phase space constraints on the scattering events, to- 
gether with the dependence of the scattering amplitudes 
on the electronic momenta, see Eqs. ^ and Q. 

The eigenstates of a translation invariant ID helical 
system can be labeled by their momenta k. If we as- 
sume that the system is TR invariant, Kramers theo- 
rem ensures that for any k there exist two degenerate or- 
thogonal eigenstates, created by the operators ip\ k and 



J +,k 

ip_ _ k , which are related by the TR operator 6, e.g. 



The kinetic energy has the gen- 



ral fori 



E 

k 
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-k)ipl_ k ip- 



(1) 



For momenta close to the Fermi momentum, k ~ kp, 
the operators k and _ k create right-moving and 
left-moving electrons, respectively, propagating with ve- 
locities ±Vp, where vp = de(k) /dk\k=k F - 

In a generic helical liquid, the electron spin component 
along a fixed z direction does not have to be a good 
quantum number. The field operators ^> a ,k of an electron 
with momentum k and spin projection a = J [, J, along the 
z-axis are related to the operators ip±,k by a momentum- 
dependent SU{2) matrix Bk, 



B h 



^+,k 



(2) 



The normalization and the orthogonality of the momen- 
tum eigenstates is reflected in the unitarity condition 
B[B k = diag(l, 1). Moreover, TR invariance entails the 
symmetry Bk = B-k which follows by comparing how 
TR affects tp a and ip a (a — ±). This also ensures that 
states created by ijy k and i/}_ a _ k at opposite momenta 

always have opposite spins, e.g., [B^B-kV = 0- As one 
consequence, elastic backscattering between such states 
is prohibited for nonmagnetic impurities. 

We take the electron-electron interaction to depend 
only on the distance between the electrons so that H lnt = 
J dxdx'U(x — x')p(x)p(x'), where p(x) = p^{x) + pi(x) 



is the total particle density. When expressed in terms of 
the eigenstates of Hq, this becomes 



H ** = \ E E U( q )[BlB k ^[Bl,B k , +q f 

kk'q a/3a'P' = ± 
X ^L#,fc-<?^< k >1pf)',k>+q, 



(3) 



where L denotes the length of the helical edge and U(q) 
is the Fourier transform of U(x). 

When considering an impurity violating translational 
invariance along the edge, we concentrate on local per- 
turbations V(x) interacting with the electron density, 
Hy = J dxV(x)p(x). Expressed in terms of ij)± ki 



H v = iJ2 E y(ki-k2)[BiB k2 r^i, kl ^M, 



feifea ap— ± 



(4) 



where V(k) is the Fourier transform of V(x). The total 
Hamiltonian H — Hq + Hi nt + By is TR invariant, as 
can be seen explicitly using the unitarity and k — > — k 
symmetry of B k . 

As we are interested in the conductance of edge 
states, we concentrate on temperatures low compared to 
the bulk gap and linearize the single-particle spectrum 
around the Fermi momentum, e(k) = vp(k — kp). We 
also make some simplifying assumptions about the form 
of B k . A /c-independent B k describes a constant rotation 
of the spin quantization axis. In this case, the spins of 
right- and left-movers are still opposite, irrespective of 
momentum. Due to the symmetry B k = B-k and uni- 
tarity, the leading terms in B k for small momenta k <C fco 
can be written as 



B k = 



1 - fc 4 /(2fc 4 ) 
P/kl 



-e/ki 
i - ky(2ki) 



(5) 



Here, fco parametrizes the scale on which the spin quan- 
tization axis rotates with k. As confirmed for certain 
microscopic models below, Eq. Km represents the generic 
form of Bk up to order (fc/fco) . Higher-order terms in 
k only give rise to subleading corrections to the conduc- 
tance. Finally, we neglect the momentum-dependence 
of the interaction and scattering potentials, and assume 
U(q) = U Q and V{k) = V . 

Combining Eqs. ([3| and ([5]), we find a TR invariant 
interaction Hamiltonian with the structure 



-Hint °C 



(fc 2 - k' 2 )U 



kk' q 



l 2 



1 P{,k+ q ^-,k'-q^+,k'^+,k 



(il> + «-> V-) + h.c. 



(6) 



H int describes backscattering of a single particle, accom- 
panied by the creation of a comoving particle-hole pair. 
The conductance correction is proportional to the rate of 
this process, which can be calculated using Fermi's golden 
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FIG. 2. Representative one-particle backscattering process to 
second order in the electron-electron interaction (wavy line) 
and impurity scattering (cross). 



FIG. 1. Dimensionless factor /(£) = f(vFkp/T) describing 
the dependence of the interaction-induced conductance cor- 
rection on the Fermi momentum, see Eq. j9|. 



rule. Of the three integrals over the momenta in the final 
state (two particles, one hole), two are canceled by energy 
and momentum conservation. At low T and kp = 0, the 
remaining momentum is of order fc ~ T/vp. The oc fc 2 
scaling of H mt then yields 5G mt °c T 5 . Interactions of 
the form (|3| also cause two-particle backscattering pro- 
cesses. However, these result in a contribution 5G oc T 9 , 
and are thus subleading with respect to the two-particle 
backscattering amplitude considered in Ref. [9]. 

The presence of impurities relaxes the requirement of 
momentum conservation in the scattering process. The 
corresponding low-energy effective Hamiltonian, applica- 
ble for kp ^ and T <C uf^f!, an d derived within 
perturbation theory in H lrA and Hy has the structure 



tion [IS] 



rreff 



(fc - k')k F U Vo 



^— ' kZvp 

kk'qq' 

+ (V^+ + h.c 



(7) 



The scattering rate now involves an additional momen- 
tum integration. Combined with the oc fc scaling of 
H^ nt , one finds SGy M oc T 4 . 

For the detailed evaluation of the conductance, the left 
and right ends (at x = ^fL/2) of the helical edge are cou- 
pled to electron reservoirs which are held at the same 
temperature T, but at slightly different chemical poten- 
tials /ii = V/2 and /ir = —V/2, respectively. In the 
clean, noninteracting limit, the conductance Go = e 2 /h 
is temperature-independent. 

We calculate the change in conductance 5G due to 
interactions and impurity scattering using perturba- 
tion theory in iJi nt and Hy. The inelastic backscat- 
tering current is determined by the transition rate 
2tt\ (/I T \i) \ 2 S(tf — ej) between initial states \i) = 

^L.feu^l^,*^ 1°) ( with ener gy e d and nnal states I/) = 

^ afl ,k fl ^a i2 ,k f2 1°) ( with ener gy e /)> weighted by ther- 
mal occupation factors. For single-particle backscatter- 
ing Ui\OLiiUf\Uj2 — — 1- The T-matrix satisfies the equa- 



T = (H int + Hy) + (H int + Hy) 



1 



£j — Ho 



f. (8) 



In the absence of interactions, Uo = 0, backscattering is 
forbidden by TR invariance. The first-order term in the 
interaction, T = H lnt , yields 



<" r , I Uq\ ( T \° r ( vp\kp\ 



5Gi nt = -rLko . . , 

n \vp ) \vpk§ 



T 



/(C) = " / dx x dx 2 (xf - x 2 ) 2 np(x 1 - Qn F {x 2 - C) 



(9) 



[1 - n F (xi +x 2 - C)][l - »f(-C)] 



Here, np(x) = l/(e x + 1) is the Fermi function, and 
the function /(C), plotted in Fig. [IJ describes the de- 
pendence on the Fermi energy. The non-monotonic be- 
havior of /(C) is a consequence of the fc-dependence of 
the matrix elements of the interaction Hamiltonian, see 
Eq. @. With the asymptotes /(£ = 0) « 306.02 and 
/(C) = (44/457r)C 6 exp(-C) for C > 1, we obtain 



5G int « 306.02 — Lk 
n 



Uo 

Up 



T 



v F k 



for kp = 0, and 
44 e 2 



407T h \Vp 



Uq \ 2 f kp \ 6 vpko 



T 



(10) 



(11) 



for «i?|fc_F| 3> T. The correction SG \ n t in Eq. (Ill is 



activated because energy and momentum conservation 
require that the counter-propagating particle in the final 
state be created at zero momentum, which is deep within 
the Fermi sea. 

In the limit uf^f! 3> T, combined processes involv- 
ing potential scattering off an impurity and scattering 
by the interaction, provide the leading contribution to 
5G. The second iteration of the T-matrix equation ^ 
yields a number of cross-terms in Hy and H lnt . A typi- 
cal contribution is shown in Fig. [2] one incoming particle 
is scattered at the impurity into a virtual intermediate 
state, and subsequently interacts with a second incom- 
ing particle. Summing all such contributions, we find the 
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effective Hamiltonian H^? int and (for Uir|feir| ^> T) 



6G Vint «1.21 x 8 4 — Ln 
h 



V U 



kp 
ko 



( T 



\v F k 



(12) 



in line with the power-counting argument given after 
Eq. Q above. In deriving Eq. (12 1 we assumed that 



the impurities are randomly positioned along the edge 
with dilute linear density n; mp , such that multiple scat- 
tering can be neglected. Specifically, interference terms 



are small for T> v F /L [16] . A comparison of Eqs. ( 11 ) 
and ( 12 ) shows that for sufficiently weak impurity scat- 
tering, (n inip /k )(Vo/vF) 2 {T/vF k ) 7 < 1CT 5 , the depen- 
dence of SG on kp at fixed T displays a minimum at some 
finite value of \kp\. 

The perturbation theory in Uo and Vq diverges if the 
initial or final state of one of the electrons is at the TR 
invariant momentum k — 0, and the intermediate state 
in Fig. [2] approaches it. This divergence is similar to the 
one in the cotunneling amplitude in Coulomb blockaded 
quantum dots, and may be treated in a similar way (see, 
e.g., Appendix C in Ref. Jl7] ) . The interactions lead to a 
finite lifetime of the intermediate particle or hole at k = 
and cut off the divergence. If that lifetime is longer than 
the time of flight L/vp, then the singularity is cut off 
by the deviation of the eigenstates from plane waves due 
to impurity scattering. In any case, higher-order terms 
in Vq and Uq regularize the divergent contributions and 
make them smaller than the results in Eqs. |9|-(12). 

So far, we have used an effective one-dimensional 
model of the helical edge. However, the edge states 
exist at the boundaries of 2D topological insulators. 
Thus, their wave functions decay exponentially on a 
momentum-dependent length scale 1/A(fc) into the bulk 
of the topological insulator, e.g., 4>k(x,y) oc e -K k )v e lk ^ _ 
As a consequence, the spin-rotation matrices involve 
convolutions over y of the two-dimensional eigenstates 
along with the 2D interaction potential U(x,y) and impu- 
rity potential V(x, y). Importantly, however, the small-fc 
behavior of Bu is always compatible with the expansion 
([5]), since the latter follows from TR invariance. 

In order to justify our effective ID model Q-Q, we 
determined the momentum-dependent rotation of the 
spin quantization axis explicitly for the BHZ model [TO] 
in the presence of Rashba SOI. This model provides 
a description of 2D topological insulators realized in 
HgTe/CdTe quantum wells. We used exact diagonaliza- 
tion to solve the Hamiltonian on a cylinder of width W, 
with periodic boundary conditions in the x direction and 
edges at y = and y = W . In the following, we shall 
outline the procedure and the results; details will be pub- 
lished elsewhere. 

The 2D topological insulators realized in HgTe/CdTe 
quantum wells can be modeled using a Hamiltonian 
H = '}2kH(k) : where H(k) has a 4 x 4 matrix struc- 
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FIG. 3. (Color online) Off-diagonal component [Bl i Bi C2 ]~ + 
of the spin-rotation matrix determined from the numerical 
solution of the BHZ model with Rashba SOI using Eq. ( 13 1. 



ture in a basis containing two spin-degenerate orbitals, 
V = {E t ,H t ,E l ,H l } [TO]. In the original BHZ model, 
H(k) is block-diagonal and does not couple spin- up and 
spin-down states. However, it was shown in Ref. |13) that 
breaking inversion symmetry by applying an out-of-plane 
electric field leads to Rashba SOI and thus a coupling be- 
tween spin-up and spin-down states. As long as Rashba 
SOI remains weak, the gapless edge states remain intact. 
For a given momentum k, H{k) has four eigenvectors 
which correspond to left-moving and right-moving states 
localized either on the upper or lower edge. These eigen- 
modes are four-component spinors in the basis V and will 
be labeled 4> a ,c,,k{x-, y), where a = ± denotes the chirality 
as before, and C = U, L labels the upper or lower edge, 
at y — W and y = 0, respectively. 

In order to translate the numerical solution of the 2D 
model into parameters of the effective ID model, we 
match the impurity scattering operator Q with a corre- 
sponding operator in the 2D system. The ID description 
of the edge state scattering can be applied if the pene- 
tration depth 1/A of the edge state into the bulk [JJ is 
small compared to the range ly of the impurity poten- 
tial. Since we are not interested in the detailed shape 
of the impurity potential, we assume that the impurity 
potential is approximately constant over the length scale 
1/A, and that ly is small compared to the Fermi wave- 
length in the edge direction. In that case, the potential 
can be treated as constant in y direction and pointlike in 
x direction. For an impurity located at position (xo,yo) 
near the upper edge, the 2D scattering operator is then 
given by H^ s — VqS(x — xo)diag(l, 1, 1, 1), where x is the 
position operator along the edge. By comparing matrix 
elements of Hb s between the eigenstates (j) a ,u.k(x, y) with 



5 



the matrix elements of Hy, we can identify 



[BlB k2 



d v4> ai ,u,ki ( x o, v)4>ai,ute {xo,y). (13) 



The result of the numerical solution is shown in Fig. [3] 
Importantly, it shows the quadratic dependence on ki 
and k,2 , in agreement with the low-momentum expansion 

In addition, we have numerically calculated the 
momentum-dependent rotation of the spin axis for vari- 
ous other 2D topological insulators with broken S z sym- 
metry, e.g., the Kane-Mele model in the presence of 
Rashba SOI, the BHZ model with bulk inversion asym- 
metry, and the model by Shitade et al. [T5] for monolayers 
of sodium iridate. We found similar results as in Fig. [3] 
in all cases. For small momenta, follows directly from 
TR invariance and translational invariance, so we expect 
it to be universally applicable for helical liquids at low en- 
ergies and weak disorder. As a consequence, the scaling 
SG oc T 4 in Eq. (12 1 can be expected to hold generically 
for helical liquids with broken S z symmetry, even if the 
spin rotation is created by other mechanisms. 
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